Algorithms for fluid-structure interaction problems arising in hemodynamics

We discuss in this thesis the numerical approximation of fluid-structure interaction (FSI) problems with a particular concern (albeit not exclusive) on hemodynamics applications. Firstly, we model the blood as an incompressible fluid and the artery wall as an elastic structure. To solve the coupled problem, we propose new semi-implicit algorithms based on inexact block-LU factorization of the linear system obtained after the space-time discretization and linearization of the FSI problem. As a result, the fluid velocity is computed separately from the coupled pressure-structure velocity system at each iteration, hence reducing the computational cost. This approach leads to two different families of methods which extend to FSI problems schemes that were previously adopted for pure fluid problems. The algorithms derived from inexact factorization methods are compared with other schemes based on two preconditioners for the FSI system. The first one is the classical Dirichlet-Neumann preconditioner, which has the advantage of modularity (i.e. it allows to reuse existing fluid and structure codes with minimum effort). Unfortunately, its performance is very poor in case of large added-mass effect, as it happens in hemodynamics. Alternatively, we consider a non-modular approach which consists in preconditioning the coupled system with a suitable diagonal scaling combined with an ILUT preconditioner. The system is then solved by a Krylov method. The drawback of this procedure is the loss of modularity. Independently of the preconditioner, the efficiency of semi-implicit algorithms is highlighted. All the methods are tested on two and three-dimensional blood-vessel systems. The algorithm combining the non-modular ILUT preconditioner with Krylov methods proved to be the fastest. However, modular and inexact factorization based methods should not be disregarded because they can considerably benefit from code parallelization, unlike the ILUT-Krylov approach. Finally, we improve the structure model by representing the vessel wall as a linear poroelastic medium. Our non-modular approach and the partitioned procedures arising from a domain decomposition viewpoint are extended to fluid-poroelastic structure interactions. Their numerical performance are analyzed and compared on simplified blood-vessel systems

modelling fluid-poroelastic media interaction

The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier-Stokes equations with the Biot system. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we first design residual-based stabilization techniques for the Biot system, motivated by the variational multiscale approach. Then, we state the monolithic Navier-Stokes/Biot system with the appropriate transmission conditions at the interface. For the solution of the coupled system, we adopt both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are considered and their convergence 1 is analyzed for a simplified problem. We compare the efficiency of all the methods on a test problem that exhibits a large added-mass effect, as it happens in hemodynamics applications.

Splitting methods based on algebraic factorization for fluid-structure interaction

We discuss in this paper the numerical approximation of fluid-structure interaction (FSI) problems dealing with strong added-mass effect. We propose new semi-implicit algorithms based on inexact block-$LU$ factorization of the linear system obtained after the space-time discretization and linearization of the FSI problem. As a result, the fluid velocity is computed separately from the coupled pressure-structure velocity system at each iteration, reducing the computational cost. We investigate explicit-implicit decomposition through algebraic splitting techniques originally designed for the FSI problem. This approach leads to two different families of methods which extend to FSI the algebraic pressure correction method and the Yosida method, two schemes that were previously adopted for pure fluid problems. Furthermore, we have considered the inexact factorization of the fluid-structure system as a preconditioner. The numerical properties of these methods have been tested on a model problem representing a blood-vessel system.&nbsp